      SUBROUTINE PrintP
        Use global
        implicit none
        integer :: i
        character(len=4)     :: NS
        character(len=20)     :: ResultFile
		
		!!! Update Pressure and speed of sound
        do i=1,nnt
		  Pres(i)=(gam-1.d0)*(U(i,3)-0.5D0*U(i,2)*U(i,2)/U(i,1))
		  V_sound(i)=Sqrt(gam*Pres(i)/U(i,1))
		enddo
        write(NS,'(i4)') istep
        NS=adjustl(NS)
        ResultFile='Step='//trim(NS)//'.txt'
        
        open(10,file=ResultFile)
          Do i=1,nnt
            write(10,100)vcor(i),Pres(i)/1.0D5,U(i,1),U(i,2),U(i,3) !,vsol(i,1),vsol(i,2)/vsol(i,1)/vcelerity(i),Var_ae(i)
          Enddo
        close(10)
        !!! x(i),p(i),rho(i),ma(i)
        100 format(f12.4,f20.4,f20.4,f20.4,f20.4)   !,f20.8,ES12.3,f10.5)
      end subroutine